On the existence of stable seasonally varying Arctic sea ice 
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Within the framework lower order thermodynamic theories for the climatic evolution of Arctic sea ice we iso- 
late the conditions required for the existence of stable seasonally-varying ice states. This is done by constructing 
a two-season model from the continuously evolving theory of Eisenman and Wettlaufer jljl and showing that the 
necessary and sufficient condition for stable seasonally-varying states resides in the relaxation of the constant 
annual average short-wave radiative forcing. This forcing is examined within the scenario of greenhouse gas 
warming, as a function of which stability conditions are discerned. 
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I. INTRODUCTION 

The recent rapid decline of the Arctic sea ice cover captures 
efforts to constrain cause, effect and sensitivity in the evolu- 
tion of the state of the system [e.g., 0,0]. A principal focus 
is the evolution of the summer sea ice minimum; if the ice 
cover vanishes in summer, some argue that this constitutes an 
irreversible change akin to a saddle-node bifurcation, whereas 
others argue that such a so-called "tipping point" is not associ- 
ated with this transition, which would then be reversible. The 
arguments involve a range of methods from theoretical treat- 
ments |[ll0-|7tl and global climate model simulations [e.g.,0- 
[lOll . to extrapolation of observations [e.g.. lTTl[T2ll . 

The heuristics of the ice-albedo feedback underlie the no- 
tion of an abrupt and irreversible transition from the perennial 
ice state to either a seasonal or ice-free sea ice state; progres- 
sive reduction of the ice cover continues due to the secular 
increase in the sensible heat of the low albedo ocean. Here we 
seek to lay bare the specific conditions under which either the 
ice-ocean albedo contrast or the ice-albedo feedback drives a 
transition to a seasonal state that is stable. 

Thorndike 101 developed an analytical model coupling sea 
ice growth to climate by calculating the annual cycle in four 
stages across a warm and a cold season. For each season he 
chose representative constant values of radiative fluxes and 
the albedos of sea ice and the ocean. His climate forcing took 
the form of a constant poleward heat transport within a range 
of which he found two stable steady state solutions and one 
unstable steady state solution. The perennial ice and ice-free 
states were stable and the seasonally-varying states, in which 
the ice vanishes each summer and returns each winter, were 
unstable. The hysteresis between these two stable states re- 
flected the presence of a saddle-node bifurcation. The origin 
of the instability was the ice-ocean albedo contrast. 

Eisenman and Wettlaufer |1|| extended the approach of 
Thorndike |0] to develop a single evolution equation for the 
state of the system forced by monthly observations of the ra- 
diative fluxes. They modeled ice-albedo feedback by treating 
the sea ice albedo as a function of sea ice thickness, transi- 
tioning continuously from that of sea ice to that of the ocean. 
In analogy with the poleward heat transport of Thorndike 
10], they forced the system with an additional heat flux IS.Fq, 
generically associated with greenhouse gas forcing, and an- 



alyzed the fixed points to their evolution equation. Upon in- 
crease of A^o they found that stable seasonally-varying states 
emerged continuously from the perennial ice states. These sta- 
ble seasonal ice states persisted for a range of AJq and then 
were lost to a saddle-node bifurcation leading to an ice free 
state. While maintaining the time dependence of the forcing, 
Eisenman Igl further simplified the model of Eisenman and 
Wettlaufer (llj to show that over a wide range of parameter 
choices the model can, among three other scenarios, capture 
the behavior shown by Thorndike |0]. Eisenman |0] provides 
a thorough summary of the models and methods used to pre- 
dict four general scenarios under which ice retreat may occur 
in a warming climate, all of which he was able to reproduce 
within the scope of his wide ranging parameter study. 

It is self-evident that the detailed construction of such con- 
ceptual models underlies the distinction between the nature 
and stability of their fixed points. The value of seeking the 
physical underpinning of this distinction is to focus think- 
ing on the essential ingredients for such approaches to cap- 
ture qualitatively distinct behaviors. That is the goal of this 
paper. In order to make the appropriate comparison with 
Thorndike |0] we construct a two season model by averag- 
ing that of Eisenman and Wettlaufer f^. In order to make 
the paper reasonably self-contained we summarize the model 
of Eisenman and Wettlaufer |0] in the next section. We then 
describe our partitioning of it into a two-season model after 
which we give the solutions in the different subseasons. The 
seasonal solution is obtained by enforcing conservation and 
energy and mass to match the ice-covered and ice-free solu- 
tions across the subintervals of time. We then find that the nec- 
essary and sufficient condition for stable seasonally-varying 
states resides in a finite time dependence of the short-wave ra- 
diative forcing. Whilst the quantitative nature of the stability 
depends on the quantitative nature of this time dependence, 
the generic structure does not. 



II. SUMMARY OF EISENMAN & WETTLAUFER (2009) 

A two season model is constructed from Eisenman and 
Wettlaufer |0] (EW09) by dividing the year into a "cold sea- 
son" and a "warm season". To insure that the reader can fol- 
low this construction we summarize EW09 here. 
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The state variable E is the energy (with units W yr) 
stored in sea ice as latent heat when the ocean is ice-covered 
or in the ocean mixed layer as sensible heat when the ocean is 
ice-free, viz.. 



E 



Ljhi 



E < [sea ice] 



CrrdHmiTmi E >0 [occan] 



(1) 



where Li the sea ice latent heat of fusion, hi its thickness, Cmi 
is the specific heat capacity of the ocean mixed layer, Hmi 
is its depth and T„ii its temperature. Ignoring salinity effects, 
the temperature T{t, E), determined by energy balance across 
the layer, is measured relative to the freezing point T,„ as 



T{t,E) = -n 



Fnit) 



hUlE - Frit) 



(2) 



where the ramp function is TZ{x > Q) ~ x and Tl{x < 0) = 
0, the thermal conductivity of the ice is ki, and the radiative 
quantities Fd and FT{t) are discussed later. 

An essential aspect of the transitions we discuss is the na- 
ture of the ice albedo feedback. Here, the Beer-Lambert law 
of exponential attenuation of radiative intensity with depth in 
a medium motivates a treatment of the dependence of the sur- 
face albedo with E using a mixture formula with a character- 
istic ice thickness ha for the extinction of shortwave radiation 



a{E) 



ami 



ami 



Ui 



tanh 



E 

Lihn 



(3) 



This describes the fraction, 1 — a{E), of the incident short- 
wave radiation Fs{t) absorbed by the ice. 

The evolution of the state of the ice (or ocean) cover is de- 
termined by the balance of radiative and sensible heat fluxes 
at the upper surface, Fd — FT{t)T{t, E), the upward heat flux 
from the ocean Fg, and the fraction of ice exported from the 
domain voTZ{—E) through a first order nonautonomous en- 
ergy balance model as 



dE 



= .f{t,E), 



(4) 



with 



fit, E)=Fd- FTit)T{t, E)+Fb + voUi-E), (5) 
where 



Fait, E) = [l- a{E)] Fs{t) - Fo{t) + AFq 



(6) 



The term Fd ~ FTit)T{t,E) is thought of as the differ- 
ence between the incoming shortwave radiation at the sur- 
face [1 — a{E)] Fs{t) and the outgoing longwave radiation 
(cx T^), augmented here by sensible and latent heat fluxes and 
an additional amount associated with greenhouse gas forcing 
AFq. The Stefan-Boltzmann equation for outgoing longwave 
radiation is linearized in the deviation of the surface temper- 
ature from the freezing point as Fo{t) + FT{t)T{t, E) which 



is known to be a reasonable approximation [e.g.,|3]. Finally, 
the seasonally varying values of Fo(t) and Frit) are deter- 
mined from an atmospheric model incorporating observations 
of Arctic cloudiness, atmospheric transport from lower lati- 
tudes and the meridional temperature gradient Finally, we 
note that the ice export vqR{—E) is typically ~ 10% yr~^, 
but the nonlinear relationship between ice thickness and ice 
growth rate highlights the possibility that in changing climates 
a time dependent value may be important in determining mul- 
tiple ice states [most recently see@, and refs. therein]. In the 
two season model we neglect ice export. 



III. A TWO SEASON MODEL FROM EW09 

The measure for the division between the cold and warm 
seasons is taken to be the downwelling shortwave radiance 
Fs{t), and we average the fluxes (all measured in Wm~^) 
involved in the seasonal evolution of sea ice over these seasons 



and 



Fs{t) 







Q<t<\ 
200 \<t<l 



104 Q<t<\ 



64 



<t<l 



Cold Season 
Warm Season, 



Cold Season 
Warm Season, 



(7) 



(8) 



wherein time is measured in years. Because Frit) does not 
change significantly over an entire season we take it as a con- 
stant 3.0 Wm^^K^^. The form of the albedo that allows one 
to study the ice-albedo feedback is given by equation (O, 
which we return to in Section lTV El but in the two season treat- 
ment when we clearly have ice we use and when we clearly 
have ocean we use ami- 

During the cold season the maximum Eq decreases even- 
tually reaching a minimum value Ei, which then increases to 
Eq again during the warm season. By imposing conservation 
of energy (and hence mass) and continuity across the seasonal 
transitions we determine the solutions E{t) for the cold and 
warm seasons for the perennial ice, ice free and seasonal ice 
states in turn. 



A. Perennial Ice Solution 

1. Cold Season 

The evolution equation during cold season is given by 

-Fo + AFo 



dE 
It 



Fn + AFn + Fu 



kiLiFnet — FbFtE 



Ft 



Ft - k,L,/E 



(9) 



(10) 
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where Fnet = ^ + + Fb, which upon defining the 
constants 



tBi'T 



gives 



A-E 



dE = Fsdi. 



B- E 

We write the imphcit form of the solution as 

Eo - E{t) 



[E{t) - Eo\ -{A- B)\n 



1 



B-E. 



(12) 



Fedt, (13) 



and because | B — Eq |:^| Eq — E{t) \ in the perennial ice 
state, then 



im - Eo] 

which we rearrange as 

E{t) ^ Eo + 
= Eo + 



A-B 



B-Eo 



FsiB-Eo) 



FBt, 



(14) 



t 



A-Eo 

kiLiFnet — FbFtEq 

hu - FtE„ 



(15) 



This expresses the stabilizing energy balance during the cold 
season; the quantity multiplying time is negative with a mag- 
nitude that increases as Eq decreases thereby expressing the 
basic physics that thin ice grows faster than thick ice during 
winter. The cold season ends at t=l/2 when E{t) reaches a 
minimum Ei given by 



El = Eq 



kiL,F, 



i-^i-^ net 



FbFtEq 



2{hL, - FtEo 



(16) 



2. Warm Season 



from which we can see that the summer sea ice cover vanishes 
when Fnet + Fnet = 0. Therefore, this condition allows us 
to determine the greenhouse forcing AFq associated with the 
vanishing of the ice cover viz.. 



AFo 



Fq + Eq 



53 Wm" 



(1-aO 



Fs 



En 



(20) 



Such a simple expression demonstrates how the imbalance of 
heat fluxes over the Arctic Ocean is compensated for by the 
growth or decay of sea ice. The condition ( l20b for the vanish- 
ing of sea ice during summer represents a turning point in this 
energy balance. 



B. Ice-Free Solution 

The other stable solution that we know exists as greenhouse 
forcing Ai^o increases is ice- free [e.g.,[llll. We construct this 
solution in the two season setting as follows. 



1. Cold Season 



The evolution equation is 



dE 

— = -rE + Fnet, 

at 



(21) 



wherer = FT/CmiHmi andagainF„et = -Fq + AFq + Fb. 
The solution is thus 



E{t) - Foe-'^* 



Fnet 



(22) 



and at the end of cold season, E(t) reaches the minimum Fi, 
which is 



Fi =Foe-^'^ + — (1-e-^'-). 

7' 



(23) 



During the warm season, the ice is ablating and hence 
r(t, F) is zero, and the evolution equation is 

dE 

— ^{l~a,)Fs-Fo + AFa+FB^Fnet, (17) 
which we integrate from t = 1 /2 to 1 giving 



Eq — El — -Fnet- 



(18) 



With the use of equation (fT6b this can be expressed in terms 
of fluxes as 



Fo 



kiLi{Fnet + Fnet) 
Fnet Ft + FbFt 



(19) 



2. Warm Season 

During the warm season the shortwave radiative flux makes 
an important contribution to the system modifying the evolu- 
tion equation as 



dE 

= -rF + F* 

dt '^^^"e 



(24) 



where F,*g( = (1 — q;„j;)F5 — Fo + AFo+Fb, and the solution 
is 



F(t) = Fie-"^* + ^(1 - e-''*). 
r 



(25) 
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For a closed energy cycle E{t) must reach Eq at the end of 
the warm season giving 



(26) 



Thus, equations i23[ and (|26] | allow us to determine Eq and 
El for the ice free state which are 



j^^^F^^^te^^l + Ke^ and 



El 



r(l + 6-5'-) 



(27) 
(28) 



Clearly the existence of ice free states requires that the mini- 
mum E{t) must be greater than 0, which here implies -Ei > 
and thus 

Fnet + Fnet^~^^ ^ which is equivalent to 

Fo-Fb- e-5-[(l _ a,ni)Fs - Fo + Eg] 



AEq > 



1 + e- 



14W m- 



(29) 



We have thus captured analytically an essential feature of the 
bifurcation diagram of EW09 and the analysis of Thorndike 
fill of a hysteresis between perennial ice and ice free states. 
These two possible stable states exist over a range of AFq de- 
termined by the minimum condition for ice free states, equa- 
tion (|29] |, and the condition determining the value at which 
perennial ice vanishes, equation (l20l i. Now we examine the 
nature of the seasonal states. 



C. Seasonal Ice Solution 

We seek to understand the conditions under which seasonal 
ice states are stable. The approach is simply to combine the 
perennial ice and ice free solutions. At the end of the warm 
season E{t) reaches the positive maximum value Eq. Before 
the end of the cold season sea ice must form such that the min- 
imum El is negative. Then, during the next warm season the 
ice ablates completely to return to Eq. Figure[T]is a schematic 
of such seasonally varying solutions. Following the notation 
above we define 

F„a = -Fn + AFn+FB, (30) 
Fr.,t = {l-a,)Fs-Fn + AFo + FB, and (31) 

= (1 - a„u)Fs -Fo + AFo + Fb, (32) 

in which the Cold Season (Fnet), the Ice Covered part of the 
warm season (Fnet), and the Ice Free Warm Season {F*^^) 
are represented. The nature of the seasonally varying solu- 
tions requires a crossover time ti from the ice free state to 
the ice covered state during the cold season, and a return to 
the ice free state at a time t2 during the warm season. Thus, 
the solution for the seasonally-varying state must be solved 




FIG. 1: Seasonally- varying solutions in the two season model. Start- 
ing at t = from a positive open ocean energy Eq at the beginning 
of the cold season, ice forms when E{t) crosses into the negative 
region at f — ti. The ice grows during winter to reach the minimum 
El at the end of the cold season, t — 1/2. During the warm season, 
E{t) becomes positive and the ice vanishes at f = t2, after which the 
ocean warms back to Eq. The net radiative forcing's Fnet , Fnet and 
F*et are described by equations l |30|32| l 




20 30 
A F„ (Wm-2) 



FIG. 2: Bifurcation diagram for the two season model as a function 
of the external, greenhouse gas forcing, heat flux AFo. Thick lines 
represent summer and winter energies calculated numerically and the 
dashed-dot lines the summer and winter energies calculated analyti- 
cally. Here, seasonal ice vanishes in a saddle-node bifurcation. 



incrementally in the four stages shown in figure [T] 
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1. 0<t<ti 



from which, with the aide of equation dJTl i. we find t2 as 



The cooUng of the ocean is governed by the evolution equa- 
tion for E{t) for the ice free state which is 



dE 
'dt 



-rE + F„ 



(33) 



We thus integrate with an integrating factor e''* from t = to 
t^ti, noting that Eit = ti) = 0, to find 



1 / r 
ti = -In 1 - —Eq 

r V Fnet 



exploiting the fact that | rEo\ |<c| Fn 



Erie 



(34) 



2. ti <t< 1/2 

Once ice has formed it grows according to the evolution 
equation we used for the perennial ice state during the cold 
season; 



dE ^ kiLiEnet — EbFtE 
~dt ~ k,Li - EtE ' 



(35) 



Integrating from t ^ ti io t = 1/2 and using A and B from 
above we find 



{El - 0) 



A-B 



B-0 



= FB{l-ti 



(36) 



and hence 



kiLiEnet/ FbFt 



hU/FT 



— Eo + -Enet- 



E ( ^ + ^° 

2 Fn.pt 



(37) 



3. l/2<t<t2 

Having passed through the minimum of E{t) the system 
begins to warm. It evolves with the albedo of sea ice and is 
described by 



, , — ^ net ; 

dt 

which we integrate from t = 1/2 to t = ^2 to find 

— El = Enet (t2 ^ 2 ) ' 



(38) 



(39) 



^ _ 1 Eq 1 Enet 
2 Enet 2 Enet 



4. t2<t<l 



(40) 



As t passes through t2 the sea ice vanishes and the radia- 
tively exposed ocean warms according to 



'4^-rE 
dt 



e: 



^-rE + {l~ami)Es-Eo + AEo+EB, (41) 
from which we find 



En 



p* 



1 



exp 



2 net 



J 7iet TP 

-r — — £/o 



1 

2 

^e: 



Eq 

Enet 
Enet 



I En 



2F, 



net 



2' F 

^ ne 



En 



e: 



(42) 



wherein we rely on the observation that 
^ (5 + 'f^ + 5^^) small throughout the entire 
range of AFq studied in EW09. This originates in the fact 
that the term describing the excess longwave radiative flux 
due to the temperature change of the ocean, —rE, is small 
relative to the net radiative flux, E*^^. However, as AFq 
becomes very large the approximation breaks down. 

The expected behavior of the system is summarized in Fig- 
ure 2. Dashed-dot lines give the behavior calculated analyt- 
ically using approach described above. Thick lines represent 
summer and winter energy values of the stable steady state so- 
lutions calculated numerically with the representative radia- 
tive flux values (the same as those in the analytical calcula- 
tion) and ice export excluded. Clearly the bifurcation points 
occur at different values of AFq, which due to the fact that 
we plot the analytical results using all of the approximations 
discussed, whereas these approximations are not made in the 
numerical calculations. However, the point is that the main 
physical interpretation is the same. We see in both cases the 
range of AFg in which the two stable solutions (the perennial 
ice and ice free states) and one unstable solution (the season- 
ally varying state) coexist. Although this bifurcation diagram 
is sufficient to interpret the stability of steady state solutions 
parametrically, we seek an expression defining the conditions 
for stability that will provide insight into the minimal physical 
constraints on the stable existence of seasonal ice. 



IV. STABILITY OF SEASONAL ICE AND THE TIME 
VARIATION OF SHORTWAVE RADIATIVE FLUX 

Having found the solutions for a two-season variant of 
EW09 we now examine their stability in the usual manner. 
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taking each in turn. 



periods and the same shortwave radiative flux, more energy is 
absorbed in the ocean than in the ice cover 



A. Perennial ice state 



First, we assess the stabiHty of the perennial ice states by 
perturbing Ei in equation ( fTSI l as E'l' ~ Ei + S' which is 
accompUshed by setting Eq to E^ + S such that \S\ <C |£'o| 
giving 



El' = Eo + S 
+ 1 



hL,F„,t~FBFTiEo + 6) 
2{hL, - Ft{Eo + S)) 

kiLiFnet — FbFtEq 



2{hL^ - FtEo) 

FTkiLi{Fnet — Fb) 



RecaUing that 

we have 

Ex =Ei+5 



1 - 



2{kiLi - FrEof 

_ kiLi{Fnet + Enet) 
EriFnet + Eb) 

EriEnet + EbY 



(43) 



(44) 



2hL,i~Fnet + Fb) 



^Ei+S', 



(45) 



and hence S' < S and the perturbation decays with time and 
that the perennial ice state solutions are stable. 



B. Ice free state 

Carrying out the same procedure with the ice-free solutions 
we find 



(46) 



Because r is positive, the perturbation decays and the ice-free 
states are stable. 



C. Seasonally varying state 

When we perturb the initial energy Eq of these states in the 
same manner we find 



(47) 



where F*^f /Fnet > 1 and thus the perturbation grows with 
time. This is just a simple but more formal demonstration 
of what we understood from independent arguments; in this 
two season model the seasonally varying solutions are unsta- 
ble. The physical origin of the inequality F*^^ > Fnet is the 
albedo difference between sea ice and water; for the same time 



D. Varying the Shortwave Radiative Flux 

The main results will now be demonstrated. The arguments 
leading to equation Wt\ show that the instability of the sea- 
sonal ice states resides in the albedo contrast between ice and 
water. The warm season heat fluxes Fn^t and E*^^, given 
by equations OTT i and ( |32] |. are the constant values during 
1/2 <t <t2 and t2 <t <1 respectively. The sole difference 
between Fnet and F*^^ resides in the albedos and hence the 
constant absorption of shortwave radiation in the ice and the 
ocean. Indeed, we emphasize that the shortwave, long wave 
and ocean heat fluxes are by prescription constants in the two 
season model. We now consider relaxing this prescription by 
allowing for some temporal variation in the shortwave radia- 
tive flux Es{t) as 



{Ene 



t2 - 1/2 Ji 



t2 



Fnet dt and 



/2 



{Enet) — 



1 - t 



2 Jt 



F* rif 
^net "'■5 



(48) 



such that we can envisage that the average warm season values 
can possibly lead to stable conditions viz.. 



(Enet ) 
{Enet ) 



< 1. 



(49) 



When perennial ice is present during the cold part of the 
seasonal ice cycle {ti < t < 1/2), the energy is governed by 
equation dSTl i as 

El = Eq + —Fnet- 

Now, as the ice warms (1/2 < t < t2) recall that the evolution 
equation is 

^ = (1 - a,)Fs{t) -Fa + AFq + Fs, (50) 

but here upon integration we maintain a time dependence in 

Esit) 



- Si = (1 - a,) / Es{t) dt + {t2 - -)F, (51) 

Jl/2 ^ 

where F = -Fq + AFq + Fb and write F{tf) = 
Jtf Esit') dt', which then leads to 



Ei^-{l-a,)T{t2)-\^t2--jE (52) 
In the same manner, upon integration of the equation for the 
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ice free system in t2 < t < 1 v/e find 

Eo = il- a„-a)[ni) - nh)] + (1 - ^2)^. (53) 

Tiius, energy conservation demands that we combine these 
two results according to equation ( |37] | to give 

^Fnet + (1 - a„u)J'{l) + F{1 - ta) 
-ia,-a^i)T{t2)+(^t2-^^F = 0, (54) 

where we maintain the ta dependence for clarity in the de- 
velopment below. First, however, we discuss the simplified 
version of the above expression 



{LW) + (1 - a„,/) J-(l) 



(55) 




t(yr) 



where (LM^) = i(F„et+^"), which is the total outgoing long- 
wave radiative flux during a year (including the ocean heat 
flux). Whence, this is a statement of the fact that the net in- 
coming shortwave radiative flux is balanced by outgoing long- 
wave radiative flux, with a contribution from the ocean. It is 
useful in many respects, one of which is that from it we can 
obtain the integrated incoming shortwave radiative flux into 
sea ice J-{t2). 

Now, we reintroduce t2 into equation ( fSSl ) in order to con- 
struct an explicit version of equation ( |49] l, the condition for 
the existence of stable seasonally varying states, which is 

^^[m-Hh)\ < j^^^ Ht.y (56) 

Upon use of the expression for J- {12), we arrive at 

, ^ 1 , {I - aMLW) + {I - a^mi)Hn 
' 2+ 2{a,-a,r.i){-{LW)) ' ^""^ 

Now we use the positivity of Fs{t), and hence the fact that 
T{t) is a monotonically increasing function of its argument, 
to deduce the requirement that 



which is equivalent to 



{LW) + {I ~ a„aW) 



(58) 



(59) 



This is a necessary and sufficient condition for the existence 
of stable seasonally varying sea ice states. In the discussion 
surrounding equation i55[ the right hand side was detailed as 
the integrated flux balance. The inequality of equation i59[ 
delivers the proviso that stable seasonal ice states exist un- 
der a restriction on the net shortwave radiance absorbed by 
sea ice during the warm season, which is clearly influenced 
by the albedo contrast. If the net radiative absorption into 
the low albedo open ocean during the warm season is not too 
large, and equation ( |49] l is obeyed, then ice can still form the 



FIG. 3: Several short-wave radiance profiles. The red curve shows 
the short-wave radiance constructed from observational data. The 
simplest treatments consistent with the model are two different con- 
stants divided at t = 0.8, given by SWOl and SW02 and SW03 has 
no variation. All have an average of 200Wm~^. 



following cold season. Over the entirety of the warm season 
(with and without ice) more net energy must be used to warm 
the ice cover and ablate it rather than be stored in the ocean, 
which imposes a time constraint. Indeed, it is evident that the 
position of t2 during the warm season is determined by this 
inequality, and we examine this in the discussion surround- 
ing figures 3 and 4. Several examples of shortwave radiances 
through the annual cycle are shown in figure 3, from which 
we determine both J-{t*) and T{t2) and compare these as a 
function of AFq in the figure 4. The green straight line shows 
J-{t2) which is to be compared with the curves of J-{t*) for 
the different shortwave radiances from figure 3. We under- 
stand then from equation ( |59] l that, for a given value of AFg, 
a stable seasonally-varying state must have J-{t*) above the 
green line describing J- [12)- Clearly, the curves associated 
with the observational shortwave radiance and the constant 
shortwave radiance (SW03) fall below the green line for this 
range of greenhouse forcing and thus describe unstable sea- 
sonal ice states. However, the curves for SWOl and SW02 
describe stable seasonal ice states from a wide range of green- 
house forcing AFq. 

Fixing all fluxes except the shortwave radiance, we inte- 
grate the model numerically and examine the Poincare sec- 
tions to study the nature and number of solutions in figure 5. 
Poincare sections are developed and analyzed as follows. We 
take one of the summer energy values En and then integrate 
for a year to get En+i- The points where En+i — En = 
represent steady state solutions in the sense that they are pe- 
riodic points of period unity. When the slope crossing zero 
is negative (positive) the solutions are stable (unstable). In 
analogy with the results of fll. for SW03, which takes a con- 
stant annual value, we find only two stable steady states, the 
perennial ice state and ice-free state, and one unstable state. 
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FIG. 4: T(ti) and J-(t*) as a function of AFq. The green curve 
shows J^(f2 ) and the other curves are the different T{t* ) 's associated 
with the short-wave radiances of figure 3. 



the seasonally-varying state. As expected from our stabil- 
ity analysis above, we find a stable seasonally-varying state 
with SWOl and SW02. Indeed, between the perennial ice and 
ice-free states, we find three seasonally-varying states, one of 
which has negative slope and hence is stable. 

Consider the unstable case for which (Fnet) < {Fnet) is 
obeyed by a perturbation to a steady seasonally-varying so- 
lution. For a perturbation leading to a lower (more negative) 
winter energy value, there will be more sea ice at the end of 
winter. Hence, there is more ice to be ablated the following 
summer and thus more survives as the system enters winter. 
The feedback is positive. Figure 6 is a schematic of the stabil- 
ity conditions and their manifestations. The slope represents 
the average total heat flux during the time period indicated 
thereby reflecting whether or not the stability condition ( |49] l 
is satisfied. Therefore, the relative size of {F*^^) and {Fnet) 
determines the stability of a solution. 



E. Sea Ice Albedo Feedback 

According to our analysis of the warm season equations, a 
time dependence of the shortwave radiance can underpin the 
stability of seasonally-varying sea ice states. Thus, it is worth 
emphasizing here that whilst an important role is played by the 
fixed albedo contrast between ice and water in the two season 
model, there is no sea ice albedo feedback. Moreover, this is 
a zeroth order model that does not include an areal fraction of 
sea ice, and hence the fixed albedo contrast cannot manifest it- 
self as a feedback. It was in this context that EW09 introduced 
the feedback formulation in equation (O. 




FIG. 5: Poincare sections of energy (W m^^ yr) for various green- 
house forcing values AFo and the short-wave radiances from figure 
3. The observed (a), two examples of simple time dependence (b), 
(c) and a seasonally constant value (d). We take one of the summer 
energy values En and then integrate for a year to get En+i- The 
steady state solutions are shown as the crossing points of the line 
En+i - En = 0. When the slope at these points is negative (positive), 
the solution is stable (unstable). 
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FIG. 6: Schematic diagram of perturbation growth (decrease) given 
to seasonally-varying solution. Both unstable (a) and stable (b) cases 
are shown. The thick lines denote steady state solutions and their 
slopes represent the average seasonal heat flux which depends on the 
presence or absence of sea ice. The dotted lines show an evolution 
of positively (E > 0) or negatively (E < 0) perturbed solutions 
starting from the beginning of the warm season t = 1/2 during a 
given year. For the unstable case (a), the perturbed solution evolves 
with the same slope as the steady state solution, but after one cycle 
the deviation of the energy from the steady state solution is positive, 
as is represented by the direction of the arrows at t = 1 /2. Therefore, 
perturbations diverge from the steady state with time. Conversely, the 
perturbed solutions in the stable case (b) converge to the steady state 
solution. 



Now, we can use equation Q and consider the ice albedo 



feedback under the assumption of a constant Fs by suitable 
modification of equation ( |48] l as 



(Fnet) 



1 



Fnet dt 



t2 — 1/2 Ji/2 

= (1 - aAv)Fs - Fo 
1 

- 1 / ^net dt 

= (1 - ami)Fs - Fo 



AFn + Fe 



AFo+Fb, (60) 



where uav is defined as 
1 



OiAV 



1 n 

t2 - 1/2 Ju2 



v{E)] dt. 



(61) 



Thus, aAv < cti, due to the ice thickness dependence 
(i.e., E dependence) of the former, but it is still larger than 
ami. Hence, we see that even with the inclusion of the 
sea ice albedo feedback, it is impossible to generate stable 
seasonally-varying states so long as the shortwave radiative 
forcing Fs is a constant as we will still have (Fnet) < {F*^i). 

However, by including both ice albedo feedback and a time 
varying shortwave flux Fs (t) we have 



(Fnet) — 



1 



t2 

-Fo 



1/2 Jl/2 ' 

AFo + Fb, 



v{E)] Fs{t)dt 



(62) 



and similar arguments as provided above upon replacement of 
ai by aAv provide 



(1 - «.) [{LW) + (1 - arni)T{l)] ^ (1 - aAv) [{LW) + (1 - 



2{a^ - ami){- (LW)) 



< 



2{aAv~ami){~ (LW)) 



(63) 



r 



Therefore, for a given AFq, the inclusion of sea ice albedo 
feedback reduces the magnitude of the stability requirement 
for seasonally-varying states. 



V. CONCLUSION 

We have addressed the underbelly of whether the vanishing 
of summer Arctic sea ice cover would indicate an irreversible 
bifurcation in the climate system. To this end, although punc- 
tiliousness was cumbersome in places, we have analyzed in 
some detail the structure of a low order two season model de- 
scribing the interaction of sea ice with the climate. The model 
was constructed by suitable simplification of the more com- 
plex single column treatment of Eisenman and Wettlaufer iQl], 
which is a non-autonomous system continuously forced by ra- 
diation climatology and other observed or inferred fluxes. The 



reason for this approach is that they found a smooth transition 
from stable perennial ice states to stable seasonally-varying 
states as greenhouse forcing (AFq) increased. Hence, they 
concluded that the loss of summer sea ice is not irreversible; 
there is no hysteresis. (Note however, the transition to an ice- 
free state is indeed a bifurcation of the saddle node type; ex- 
hibiting substantial hysteresis with the perennial ice state). On 
the other hand, Thorndike jj] considered only a two season 
theory forced by constant values and found no stable seasonal 
ice. Therefore, we sought to understand the minimal physical 
conditions for the existence of stable seasonally varying ice 
states and focused our search on the basic difference between 
two season and continuously forced theories that could lead 
to robust qualitative distinctions in behavior. Such analysis is 
not possible with more complex models. 

First, the approach allowed us to capture analytically an es- 
sential feature of the bifurcation diagram of EW09 and the 
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analysis of Thorndike |0], both of which showed a hystere- 
sis between perennial ice and ice free states. We found that 
these two possible stable states exist over a range of green- 
house forcing. Simple analytic expressions were found for the 
value of AFq beyond which summer ice would vanish (equa- 
tion l20b . and below which ice free conditions can no longer 
persist (equation|29]l. 

Second, we examined the nature of the seasonal ice states 
and determined that they are unstable in this two season model 
the moment the summer ice vanishes. Equation ( |47l ) describes 
the imbalance driving the instability as residing in the albedo 
difference between sea ice and ocean which, for the same time 
periods and the same shortwave radiative flux, is responsible 
for more energy being absorbed in the ocean than in the ice 
cover Thus, the subsequent sensible heat of the ocean is too 
large to allow ice to form the following winter Whilst we 
have not examined all physically realistic two season models 
it does not appear possible to stabilize season ice with constant 
forcing, as shown in figures 3 and 4. 

Third, the form of equation ( l47b suggests that the simplest 
means of stabilizing seasonal ice states is to introduce a sea- 
sonality of the shortwave radiance Fs{t). This led to the 
ansatz (equation |49]l suggesting that during the warm season 
reversing the ratio of the radiative balance between the ice 
covered and ice free states under transient forcing can have 
a stabilizing influence on seasonal ice states. Pursuing this 
ansatz lead to an expression (equation |59]l describing the nec- 
essary and sufficient flux conditions for such stable states. 
Simple examples are given in figures 3 and 4, demonstrated 
numerically in figure 5 and shown schematically in figure 6. 

Fourth, two interesting results were found upon inclusion of 
the sea ice albedo feedback as characterized by equation Q. 
We determined that the stability condition ( |59] l is made less 
stringent by the implicit seasonal dependence of the albedo 



of the ice-ocean system, as shown in equation ( |63] |. More- 
over, despite the presence of the sea ice albedo feedback, it is 
impossible to generate stable seasonally-varying states if the 
shortwave radiative forcing Fs is a constant. This codifies 
the generally understood importance in the real system of the 
timing of ice advance and retreat relative to solar insulation. 

Finally, it is of broad interest to understand the nature of ice 
decay in the Arctic. It is extremely difficult to use comprehen- 
sive GCMs to understand the qualitative distinction between 
the range of scenarios proposed [see|6l and refs. therein], and 
extrapolation of observations on seasonal time scales is un- 
wise lfl3ll . The simplest theoretical approaches in this field 
began with the two season model of Thorndike which 
predicts that once summer ice vanishes there is an irreversible 
change to the ice free state. Whereas, EW09 and variants of 
continuously forced simple models suggest that such a change 
is reversible. Having found the reason for this distinction, 
residing in the time variation of shortwave radiative forcing 
during the warm season, our analysis provides a threadbare 
framework to examine leading order effects on the nature and 
number of qualitative transitions in the state of the ice cover. 
As such it is complimentary to the simplified version of EW09 
studied by Eisenman and, in the spirit of Thorndike 14], it 
provides a number of simple expressions trivial to use in the 
examination of the sensitivities of qualitative transitions. 
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